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Abstract 

In any geometrically nonlinear, isotropic and quadratic Cosserat micropolar extended contin¬ 
uum model formulated in the deformation gradient field C := V(/3 : 0 —> GL~'’(n) and the 
microrotation field i? : 0 —)■ S0(n), the shear-stretch energy is necessarily of the form 

:= fi ||sym(i?^F — 1)||^ -|- /Xc ||skew(i?^F — 1)||^ . 

We aim at the derivation of closed form expressions for the minimizers of in 

SO(3), i.e., for the set of optimal Cosserat microrotations in dimension n = 3, as a function of 
F G GL''’(n). In a previous contribution (Part I), we have first shown that, for all n> 2, the 
full range of weights /r > 0 and yic > 0 can be reduced to either a classical or a non-classical 
limit case. We have then derived the associated closed form expressions for the optimal planar 
rotations in SO(2) and proved their global optimality. In the present contribution (Part II), 
we characterize the non-classical optimal rotations in dimension n = 3. After a lift of the min¬ 
imization problem to the unit quaternions, the Euler-Lagrange equations can be symbolically 
solved by the computer algebra system Mathematica. Among the symbolic expressions for 
the critical points, we single out two candidates rpolar)^ (F) G SO(3) which we analyze and 
for which we can computationally validate their global optimality by Monte Carlo statistical 
sampling of SO(3). Geometrically, our proposed optimal Cosserat rotations rpolar)^^^^(F) act 
in the “plane of maximal strain” and our previously obtained explicit formulae for planar 
optimal Cosserat rotations in SO(2) reveal themselves as a simple special case. Further, we 
derive the associated reduced energy levels of the Cosserat shear-stretch energy and criteria 
for the existence of non-classical optimal rotations. 
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1 Introduction 


In this second part (Part II) of a series, we consider the weighted optimality problem for the 
Cosserat shear-stretch energy : SO(n) x GL''"(n) —)■ K(}", 

:= ^ ||sym(i?^F - 1)||^ + /Xc ||skew(_R^F - 1)||^ . (l.I) 

The arguments are the deformation gradient field F :=S/ip : il. ^ GL~'’(n) and the microrotation 
field i? : n —>■ SO(n) evaluated at a given point of the domain fl. This energy arises in any 
geometrically nonlinear, isotropic and quadratic Gosserat micropolar continuum model. Note that 
it is always possible to express the local energy contribution in a Cosserat model as IT = W{U), 
where U := is the first Cosserat deformation tensor. This reduction follows from objectivity 
requirements and has already been observed by the Cosserat brothers [171 p. 123, eq. (43)], see 
also [28) and |52j . Since U is in general non-symmetric, the most general isotropic and quadratic 
local energy contribution which is zero at the reference state is given by 

. ( 1 . 2 ) 

“volumetric energy” 

The last term will be discarded in the following, since it couples the rotational and volumetric 
response, a feature not present in the well-known isotropic linear Cosserat models[^ 

Let us now proceed to the primary objective of our present contribution 

Problem 1.1 (Weighted optimality in dimension n = 3). Let p > 0 and pc > 0. Compute the set 
of optimal rotations 

argmin W^_^^(i?;F) := argmin |/r ||sym(i?^F — 1)||^ -I- ||skew(i?^F — 1)||^| (1.3) 

flG SO(3) ’ i?G SO(3) ^ 

for given parameter F € GL“''(3) with distinct singular values cti > (T 2 > us > 0. 


11 _ 11 2 II _ 11 2 

/X sym(t/ — 1) -b skew(C/— 1) 


“shear—stretch energy” 


We use the notation sym(X) := + X'^), skew(X) := ^{X — X'^), (X, Y) := tr and 

we denote the induced Frobenius matrix norm by ||X||^ := (^X, X) = X)i<i j<n ^ij- mechanics 
applications, the weights p > 0 and /Xj, > 0 can be identified with the Lame shear modulus /x > 0 
from linear elasticity and the so-called Cosserat couple modulus pc > 0. The parameter A in the 
most general form of the energy (1.2) can further be identified with the second Lame parameter. 
Note that the interpretation of the Cosserat couple modulus Pc is somewhat delicate, see, e.g., |5D] . 
which is one of the fundamental motivations for this second contribution in a series. 


In Part I of this paper [53], we have proved a still surprising reduction lemma [531 Lem. 2.2, p. 4] 
for the material parameters (weights) p and pc which is valid for all space dimensions n > 2. This 
lemma singles out a classical parameter range Pc ^ p > 0 and a non-classical parameter range 
Pc > p > 0 ior p and pc and reduces both ranges to an associated limit case. The classical limit 
case is given by (/x, Pc) = (1,1) and the non-classical limit case is given by (/x, Pc) = (1,0). We then 
apply the parameter reduction [^ Lem. 2.2, p. 4] to Problem |1.1| and solve it in dimension n = 2. 
This allows us to discuss the optimal planar Cosserat rotations and we observe that the classical 
and the non-classical parameter ranges for p and pc characterize two substantially different types 
of optimal Cosserat rotations. 

To explain this difference, we first need to introduce the polar factor Rp{F) G SO(n) which 
is obtained from the right polar decomposition F = Rp(F)U(F) of the deformation gradient 
F G GL~^(n). Here, U(F) := 'J F'^F G PSym(n) denotes the positive definite symmetric right 
Biot-stretch tensor. We recall that the eigenvalues of [/ G PSym(ri) are by definition the singular 
values (Ti > (72 > ca > 0 of the deformation gradient F G GL~^(n). 


In the classical parameter range Pc > p > 0, the polar factor Rp admits a variational characteriza¬ 
tion which is noteworthy in its own right: namely, for all n > 2, it is the unique minimizer for (l.I) 
as a generalized version of Grioli’s theorem shows, see [531^133], or (53] Cor. 2.4, p. 5]. This 
variational characterization of the polar factor inspired us to introduce the following 


'^The Cosserat brothers never proposed any specific expression for the local energy W = W(U). The chosen 
quadratic ansatz for W = W(U) is motivated by a direct extension of the quadratic energy i n the linear theory of 
Cosserat models, see, e.g. 1431 Ibln l68l . We consider a true volumetric-isochoric split in Section : 


3.4 
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Definition 1.2 (Relaxed polar factor(s)). Let /i > 0 and fj,c > 0. We denote the set-valued 
mapping that assigns to a given parameter F G GL^(n) its associated set of energy-minimizing 
rotations by 

rpolar^,^^(J^) := argmin . 

fle SO(n) 


In dimensions fc = 2,3, we denote the associated optimal Cosserat rotation angles by C 

(—7r,7r]. More generally, in what follows, we shall denote the rotation angle of the (absolute) 
rotation field R G SO(fc) by a G (—7r,7r] and the rotation axis by r G By C K”, we 

denote the unit n— 1-sphere. In dimension k = 3, we use the well-known axis-angle parametrization 
of rotations which we write as [a, r^]. 

Since the classical parameter domain /ic > /i > 0 is very well understood by now, we can focus 
on the non-classical parameter range > fj,c > 0 m our efforts to solve Problem 0 Here, the 
parameter reduction lemma [291 Lem. 2.2, p. 4] allows us to restrict our attention to the non- 
classical limit case (/r,^c) = (1,0), because it shows the equivalence 

argmin P’) = arg min IPi_o(l?; (1-4) 

fiGSO(n) ReSO(n) 

for all n > 2. On the right hand side appears the rescaled deformation gradient F^^^^ := ■ F G 

GL~''(n) which is obtained from F G GL~'’(n) by multiplication with the inverse of the induced 
sealing parameter '■= > 0. We note that we use the previous notation throughout the 

text and further introduce the singular radius On := . 

It follows that the set of optimal Gosserat rotations can be described by 

rpolar^,^,(l^) = rpolari (1.5) 

for the entire non-classical parameter range /r > /rj, > 0. This simplifies our main objective Problem 
0 considerably, since it suffices now to solve 

Problem 1.3 (Weighted optimality in the non-classical limit case (/r, pic) = (1, 0)). Let /i > 0 and 
pic > 0- Compute the set of optimal rotations 

argmin Wi.o(i?; f) := argmin ||sym(i?'^F — 1)||^ (1-6) 

flGSO(3) i?GSO(3) 

for given parameter F G GL“''(3) with distinct singular values cti > <72 > cr 3 > 0. 


Regarding our Problem 1.3 at hand, we will see in Section that there are in general two energy¬ 
minimizing solutions with a certain symmetry. They both have the same rotation axis but differ by 
the sign of their respective rotation angles which allows us to select the corresponding branches by 
that sign. Accordingly, we introduce the notations rpolar)^_^^(F) and a^^^^{F). Loosely spoken, 
we will see that the optimal Gosserat rotations coincide with the polar factor Rp in a certain 
compressive regime for F G GL^(3), but deviate in a certain expansive regime. We shall precisely 
characterize this in terms of the singular radius Such a material behavior is commonly 

referred to as a tension-compression asymmetry which is an interesting natural phenomenon studied 
in the material sciences, see, e.g., |3^[33] and |5T] for experimental studies of nickel titanium (NiTi) 
shape memory single crystals for a glimpse on this broad subject]^ 


Problem 1.3 is a minimization problem on the matrix Lie group SO(3) of special orthogonal matrices 
parameterized by the deformation gradient F G GL’'’(3) in the identity component of the general 
linear group. Although it is not our duty, we want to point to some valuable introductory references 
to this subject. An excellent general reference for minimization problems on manifolds is the text 
by Absil, Mahony and Sepulchre [T]. There, also numerical solution approaches are presented. For 
an introduction to Lie groups and matrix groups, we refer to, e.g., [SUIT] and [H]. For compact 
Lie groups and their representation theory, see, e.g., |42] and |T0]. There is also a growing body 
of closely related work treating minimization problems on matrix groups and Grioli’s theorem in 
a similar context 


^We do not claim that such materials can actually be realistically modelled as a Cosserat continuum, although 
it is not impossible. 
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Instead of turning towards the solution of Problem 1.3 right away, we take a step back and notice 
that there is still another opportunity for simplification which reduces the space of parameters 
F e GL''"(3) to the space of ordered singular values CTi > 0-2 > 0-3 > 0 of F. This can be achieved 
by a principal axis transformation which introduces a relative rotation R and allows us to introduce 


Definition 1.4 (Cosserat shear-stretch energy for the relative rotation R). Let /i > 0, /Tc > 0 and 
let D := diag(tTi, (72, ( 73 ) with cti > <72 > era > 0 the singular values of F G GL“''(3). The energy of 
the relative rotation R G SO(3) is given by 


:= W^^^SR^-D) := yi sym(i?D - 1 ) 


skew(i?Zl — 1 ) 


(1.7) 


This transformation is described in Section |2] and leads us to the reduced 


Problem 1.5 (Optimality of relative rotations in dimension n = 3). Let /i = 1 and fic = 0. 
Compute the set of optimal relative rotations 


argmin Wip{R-, D) 

ReSO(3) 


arg min 
flG SO(3) 


sym(i?diag((7i,t72,(73) 



( 1 . 8 ) 


for a given diagonal matrix D diag(cri, t 72 , < 73 ) with cti > (72 > (73 > 0 the ordered singular values 
of the deformation gradient F G GL'*’(3). 


In this text, we strive to mark quantities related to relative rotations with a “hat”-symbol, e.g., 
we write R G SO(3). Further, we note that although, for now, we explicitly exclude the case 
of multiple singular values of F from our analysis, there is no major obstruction. The technical 
treatment would, however, clutter our exposition of the basic mechanisms which we want to distill 
here. 

At present, an explicit formal solution for the three-dimensional problem (let alone the n- 
dimensional problem) seems out of reach for us. We have, however, successfully computed explicit 
formulae for the critical points of the Gosserat shear-stretch energy by the use of computer al¬ 
gebra from which we have determined optimal solutions. For this approach to succeed, we first 
lift the Cosserat shear-stretch energy expressed in principal axis coordinates to the sphere of unit 
quaternions C El and subsequently apply the Lagrange multiplier technique for minimization 
with equality constraints, see, e.g., [39]. The unit quaternions form a two-sheeted cover of SO(3) 
and allow for a convenient representation of rotations in three-space. For a preceding successful 
application of quaternions to represent the rotational degress of freedom in Cosserat theory, see, 
e.g., 1^. A highly interesting recent approach to Cosserat shell theory which also uses quaternions 
is based on geodesic finite elements, see m and [ssiEnj. 

This paper is now structured as follows: in Section we introduce the lift of the Cosserat shear- 
stretch energy from SO(3) to the sphere of unit quaternions C El = K^. We then state the 
corresponding Euler-Lagrange equations and present the energy-minimizing solutions. The com¬ 
plete set of critical points computed by Mathematica |5S] is provided in Appendix [A| In Section]^ 
we present a geometric interpretation of the optimal Cosserat rotations rpolar)^^^(F') in terms of 
the maximal mean planar stretch for the entire non-classical parameter range ^ > /ic > 0 . 

This leads us to introduce a classical and a non-classical domain for the parameter F G GL~'’(3) for 
which we also derive some interesting alternative criteria. This illuminates the bifurcation behav¬ 
ior of rpolar^^^^(E). Further, we compute the associated reduced energy levels kF™))^(E) for the 
Cosserat shear-stretch energy. Then in Section [^ we shed light on our methodology for the anal¬ 
ysis of the critical points and the experimental computational validation of the energy-minimizing 
Cosserat rotations using statistical (Monte Carlo) methods. Finally, we summarize our findings in 
a short conclusion presented in Section 


2 Solvable Euler-Lagrange equations: transformation, lift 
and Lagrange multipliers 

In this section, we use a classical result from the representation theory of compact Lie groups 
to cast the reduced minimization problem stated as Problem |1.5| into a form which allows us to 
symbolically compute explicit expressions for the critical points using Mathematica. 
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It is well-known that the Lie group of unit quaternions C El = is closely related to the matrix 
group of rotations SO(3), see, e.g., [49] or [SH Chap. 9]. More precisely, the unit quaternions 
form a double cover of the matrix group SO(3). For a general introduction to analysis on smooth 
manifolds which includes smooth coverings, see, e.g., |47j . For a dynamical systems approach to 
quaternions, see, e.g., [72] which nicely demonstrates the usefulness of quaternions for mechanics 
applications with constraints. A more algebraic approach to quaternions with historical remarks 
is given in [25], and, finally, for those who enjoy the classics, see [37] and [38] . 

2.1 Transformation into principal directions 

In order to reduce the parameter space GL'''(3), we use the (unique) polar decompositiorj^ 
F = Rp [/ and the (non-unique) spectral decomposition of C/ = Vf'^F S PSym(3) given by 
U = QDQ^, Q G SO (3), and expand 

R^F = R^RpU = R^RpQDQ^ . (2.1) 

Here, the diagonal matrix D = diag((Ti, ( 72 , 0 - 3 ) contains the eigenvalues of U on its diagonal. These 
are precisely the singular values of F G GL“''(3). In fact, this is a particular form of the singular 
value decomposition (SVD). If F has only simple singular values, then it is always possible to 
choose the rotation Q such that an ordering cti > (T 2 > ua > 0 is achieved. 

Exploiting that Q G SO(3), it is now possible to carry out a transformation of the Gosserat shear- 
stretch energy into principal axis coordinates - essentially due to isotropy of the energy. For the 
actual computation, note first that 

(5'^(sym(i?^F) - 1)Q = Q'^ (sym(i?'^ Rp QDQ^) - 1)) Q 

= sym((5^i?^ Rp QDQ^Q — Q^Q) = sym((5^i?^ Rp QF — 1) = sym(i?F — 1) . (2.2) 

=: R 


In the process, it is natural to introduce the rotation 

R := Rp Q (2.3) 

which acts relative to the polar factor Rp in the coordinate system induced by the columns of Q, 
i.e., in a positively oriented frame of principal directions of U. This interpretation is also nicely 
illustrated by the inverse formula 

R = (qRQ'^ Rp^) ^ = Rp QIFq'^ (2.4) 


step is to insert the transformed symmetric part ( 2 . 2 ) into the definition of 


which allows to recover the original absolute rotation R from the relative rotation R. Our next 

(2.5) 


;-F) = sym(i?^ F — 1 ) = sym(i?^ F — 1)(3 = sym(i?F — 1 ) 


where we have used that the conjugation by preserves the Frobenius matrix norm. 

This is a promising simplification of the Gosserat shear-stretch energy, because it reduces the 
dimension of the parameter space from dimGL~*’(3) = 9 to only 3 parameters. However, we still 
have to account for the non-uniqueness of Q. To this end, we introduce the following symmetric 
rotation matrices 


Qi := 1 , Q 2 := diag(l,-l,-l), Qa := diag(-l, 1 ,- 1 ), (54 := diag(-l,- 1 , 1 ) , 

and collect them in a set S ;= {Qi, < 52 , Qs, Q4} C SO(3). This set forms a discrete subgroup 
of SO(3) which is isomorphic to the Klein four-group F 4 = 1? x Z^, as is easily inferred by a 
comparison of the multiplication tables. 

^For an introduction to the polar and singular value decomposition, see, e.g., m and for recent related results 
on variational characterizations of the polar factor Rp(F), see |44l 1691 ItT) and references therein. 
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Remark 2.1 (Uniqueness of the factor Q). Let ai > a 2 > cr^ > 0 be the ordered eigenvalues of 
U := VF'^F and let D := dietg{ai,a 2 , < 73 ) ■ It is well-known that the factor Q G SO(3) in the spectral 
decomposition U = QDQ^ is only determined up to the choice of a right handed orientation of the 
uniquely determined orthogonal eigenspaces of U. This corresponds to the products QS, S G S, 
which represent all of these possibilities. 


It is easy to see that for any possible choice of right handed orientation encoded by S' G 5, we 
obtain the same energy level 


sym(RD — 1) 


2 


sym{RSDS'^ — 1) 


sym{S'^RSD - 1) 


which implies 

argmin ;-D) = {S^RS | i? G argmin IUi_o(-R^ ; ^) and S G S} . 

flGSO(3) flGSO(3) 


( 2 . 6 ) 


(2.7) 


Thus, S is a symmetry group of the set of energy-minimizing rotations which acts by conjuga¬ 
tion. The previous analysis reveals that the non-uniqueness of Q G SO(3) is not an issue for the 
minimization problem, since all possible choices QS, S G S, lead to the same energy levelj^ 

Without any loss of generality, we may henceforward focus on the solution of 


arg min 
hgSO(3) 


syui{RD — 1) 


2 


argminIUi_o(^'^ ;^) • 

ReSO{3) 


( 2 . 8 ) 


This proves the reduction of Problem o to the minimization problem described in Problem |1.5| 
in Sectionj^for the non-classical limit case (/i, Hc) = (1,0). The same principal axes transformation 
can also be carried out for arbitrary values of p, and pc which gives rise to Definition |1.4| 

In what follows, we denote the rotation angle of the (absolute) microrotation field R G SO (3) by 
a G (—7r,7r] and the axis of the rotation by r G S^, where S” C denotes the unit n-sphere. 

This leads us to the axis-angle representation of a rotation which we write as \a, r'^]. In what 
follows, we work with different parametrizations of the group of rotations SO(3) simultaneously. 
Thus, we introduce the symbol = in order to identify rotations in SO(3) which are described with 
respect to different parametrizations of SO(3). For example, we might write for the relative rotation 
R = [/3, (fi, f 2 , ra)] and for a unit quaternion q G describing R G SO(3), we have q = R = —q. 
We see that, in general, this binary relation is not unique since the parametrizations need not be 
one-to-one. 

The symmetry group S := {Qi, Q 2 , Q 3 , Q 4 } hints at the structure of the set of optimal relative 
Cosserat rotations. In our previously introduced notation, we find: 


QiRQi 

QiRQs 


= 

/ 3 , (D,r2,r3) 


Q 2 RQ 2 = 

/ 3 , {fl,-f2,-f3) 

= 

/ 3 , (-^ 1 ,^ 2 ,- 

fa) 

, QJRQa = 

/3, (-fi,-r2,f3) 


(2.9) 


We observe that for rotations about the coordinate axes, i.e., with f = e^, n = 1,2, 3, the rotation 
axis f is either left invariant or negated. The latter is equivalent to the negation of the rotation 
angle fi. 


2.2 Lifting the minimization problem to 

The unit quaternions can be identified with the three-sphere := {g G El | |g| = 1} which we shall 
consider as a submanifold of the ambient coefficient space of the quaternion division ring H. Let 
us choose the coordinates {w, x, y, z) G K^, i.e., we write quaternions as q = w ix jy kz G M.. 

In order to cast the minimization problem into a form which lends itself to the derivation of a closed 
form solution, it is helpful to simplify the domain of minimization, i.e., to choose a well-adapted 
system of coordinates. We achieve this by lifting the Cosserat shear-stretch energy from SO(3) to 

consistent choice of Q{F) S SO(3) for different values of G GL+(3) is certainly to be advised for the 
numerical computation of a field of minimizers rpolarjj;^^^ (ir(a;)) depending on a; S £7. The inversion formula l |2.4[ l 
explicitly depends on the choice of Q and is sensitive to flips of the subspace orientation Q 1 —^ QS, S G S. 
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the covering space given by the sphere of unit quaternions C The principal idea is then to 
extend the covering map from to the ambient space and to apply the Lagrange multiplier 
rule with the constraint function g{q) := |g|^ — 1 = 0. This approach leads to minimizers in the 
submanifold of unit quaternions q G which project to energy-minimizing rotations under the 
well-known covering homomorphism 

/1 — 2{y'^ + z'^) 2{xy — wz) 2{xz + wy) \ 

TT : SO(3), T^iq) = I 2{xy + wz) 1 — 2{x‘^ + 2{yz — wx) j . (2.10) 

y 2{xz — wy) 2{yz + wx) 1 — 2{x‘^ + y"^)J 

In order to make our procedure explicit, let us first consider the case of arbitrary smooth energies 
W : SO(3) ^ K. 

Lemma 2.2. Any smooth energy W : SO(3) —>■ K admits a lift to a smooth energy W'^ : —>■ K 

S3 

SO(3) 

such that minimizers of are projected to minimizers ofW, i.e., 

7r( argminlT**((7) ) = argmin tT(i?) . (2.11) 

gGS3 flGSO(3) 

Proof. The covering map tt : S^ —> SO(3) defines a surjective Lie group homomorphism with 
kerTT = {1,-1}, see, e.g., |34j . This implies that the unit quaternions form a two-fold cover of 
SO(3). In particular, the Lie group homomorphism tt is a local diffeomorphism when restricted 
to a sheet of the covering and maps critical unit quaternions in S^ to critical rotations in SO(3). 
By definition W'^{q ; F) : S^ x GL''‘(3) —>■ Kq , := IT o tt is a lift of the Cosserat shear-stretch 

energy to the covering space S^. Smoothness of IT® is obvious since the composition of smooth 
maps is smooth. ■ 



For any R G SO(3) there exists a g G which represents this rotation as R = 7r(g) G SO(3). 
However, this representation is only unique up to antipodal identification, i.e., q and —q represent 
the same rotation: 7r(g) = R = 7r(—g). We further note that tt can be symbolically evaluated for 
all g G El which induces an extension. 

As previously defined, the covering map tt is only defined for unit quaternions g G S^. However, in 
order to apply the Lagrange multiplier theorem we have to extend it to a suitable neighborhood 
in the ambient space. To this end, we introduce the punctured space of non-zero quaternions by 
H := EI\{0} = K"^\{0}. Then, identifying ^ ]j9 concatenation of rows, allows us to 

consider tt : H — >■ map tt : El —)■ which leads us to the following matrix representation 

of the derivative 


^iw,x,y,z)’^{q{w,x,y,z)) 



-2z 

22/ 

2z 

0 

—2x 

-22/ 

2x 


T 

0 

21/ 

2z 

22/ 

—4x 

—2w 

2z 

2w 

—4x 


-4y 

2x 

2w 

2x 

0 

2z 

—2w 

2z 

-4y 


\-4z 

—2w 

2x 

2w 

-4z 

22/ 

2x 

22/ 

oj 



It is not hard to infer that Dq7r(g) is of rank 4 for all g G BI. Hence, the implicit function theorem 
ensures that tt : BI —>■ K^xs jg g ^Q^gj diffeomorphism from the punctured ambient space BI of the 
unit sphere to its image 7r(BI) C 


Definition 2.3 (Extension of the lifted energy). The extension of the Lie group homomorphism 
TT : S3 — >• SO(3) to BI given by tt :W ^ M^xs an extension of the lifted energy to the 

ambient space Bl 

irforBI^K, lTfo(g;D) := ||sym(7r(g)D - l)f . (2.13) 
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Let us abbreviate R{q) ■= Trjga (q). It is precisely the restriction of the lifted energy to the unit 
quaternions for which the Cosserat shear-stretch energy of the relative rotation is well-defined 


^1,0 ( 9 ;^) 


Wi^oiR{q)\D) . 


This extension is simply a mathematical construction, i.e., for g G H \ the lifted energy 
loses its original interpretation as a shear-stretch energy. Further, we note that the choice of 
extension is not unique, but the solutions to the Euler-Lagrange equations do not depend on the 
particular extension|^ 

Definition 2.4 (Lagrange function). Consider the constraint function g : El —>■ K, g{q) := \q\^ — 1. 
The Lagrange function for g : H —)■ M zs given by 

Li^o : H X K K, £ 1 , 0 ( 9 , A ;D) := wfQ{q;D) - Xg{q) . 


Clearly, g{q) = 0 if and only if 9 G C H which leads us to our final reformulation of the 
original Problem o in terms of quaternions describing relative rotations, namely 

Problem 2.5 (Lagrange multiplier formulation). Compute the critical points of the Lagrange 
function 

£ 1 . 0 ( 9 , A ;£)) = ||sym( 7 r( 9 )D - l)f - A (| 9 |^ - 1 ) (2.14) 

and determine the energy-minimizing solutions. 


The Lagrange function is polynomial. Thus, the application of the Lagrange multiplier technique 
leads to an algebraic problem for the Euler-Lagrange equations which we investigate next. 


2.3 Euler-Lagrange equations, critical points and optimal solutions 


In what follows, a shorthand notation is helpful, so let us introduce 

Sij 1 = (Tj (Tj and d^j (Jj — 9 — t, ^ * 


Towards a derivation of the Euler-Lagrange equations in quaternion representation, we first com¬ 
pute the product 


/I - 2 (7/2 + z^) 

2{xy — wz) 

2(xz -b wy) \ 


0 

0\ 

7 r(q(w,x,y,z)) B = 1 2 (xy-h wz) 

1 — 2 (a;2 -b z2) 

2{yz — wx) 

0 

CT2 

0 

y 2(xz — wy) 

2{yz -b wx) 

1 - 2 (a;2 -b j/^)/ 

VO 

0 

CT 3 / 


/(Ti (1 - 2 -I-z^)) 2a2{xy-wz) 2a3{xz-\-wy) \ 

I 2 ai(xywz) (J2 (l — 2 -I-z^)) 2a3{yz — wx) 

\ 2 ai{xz — wy) 2a2{yz-\-wx) (T3 (l — 2 (a;^-I- 


From this, we infer the symmetric part 


sym {TT{q)D) 


CTi (1 - 2(i/2 + z^)) 

S12 xy -I- di2 wz 
S31 xz dsi wy 


S 12 xy di 2 wz 
02 (1 - 2{x'^ -j- z^)) 
S23 yz d23 wx 


S31 a;z -I- dsi wy \ 
S23 yz d 23 wx 

0-3 (1 - 2 (x^ +y^))J 


Observing that sym {TT{q)D — 1 ) = sym ( 71 ( 9 ) 1 ?) — 1, we can compute the square of the Frobenius 
norm. This yields the following explicit expression for the Lagrange function Li,o : H x M —)■ K: 


£1,0(9, A ;D) = ((Ti(1-2j/2_2z 2) - l)" + {a2{l-2x^-2z^) - l)" + {a 3 {l- 2 x^- 2 y^) - l)" 
-I- 2 (^{si2xy di2wzf‘ + {S31XZ + dsiwyf {s 23 yz -b d 23 wxf‘^ 

- \{w'^ + x^ + y^ + z^ - 1 ) . 


^Alternatively, one may use, e.g., the following extension which yields pairwise orthogonal columns 

' + — ip' — 2{xy — wz) 2(xz + wy) 

7r'(g) := I 2(xy + wz) w^ — x^-h y^ — z^ 2(yz — wx) 


2(xz — wy) 2 (y 2 : + wx) w^ — x^ — y^ -h z'^ 


The restrictions Trjgs = Tr'jga to the sphere of unit quaternions are identical. 




Let D — diag(cri, (J 2 , ua) be given. Then a critical tuple of coefficients {w, x, y, z, A) for the Lagrange 
function Li^q satisfies the Euler-Lagrange equations in quaternion representation, i.e., 

^1,0 Z/i A , D') 0 . (2.15) 

After a lengthy computation in components (for which we have used Mathematic a), one obtains 
an explicit form of the Euler-Lagrange equations for Li q which is equivalent to the following 
parameter-dependent system of polynomials 

n ( j2 2 , .2 2 , .2 2 A A 

0 = w • ( ^23 X -I- dai y + ai2 z ~ ^ j 

0 = X • ^^23 -|- 4(a2 -|- (73) + (4(T3 -|- S12) + (4(72 + S31) — ((^23 -|- (S23 — 2)s23)) — 

0 = y • -I- 4 (ai -I- (7?) -|- (4(7? -|- sia) -|- (4(7? -|- sfa) x^ - {dh + (ssi - 2)s3i)) - 

0 = Z ■ ^di 2 W -|- 4((7i -|- ( 72 ) Z -|- ( 4(72 + S 31 ) X -|- (4(7x -|- S 23 ) y — (di 2 + (si 2 — 2)si2)) — —^ 

0 = -I- x^ -I- y^ -I- z^ — 1 . 


In general, solution sets of polynomial systems over the field of complex numbers C define complex 
varieties which intuitively can be regarded as almost-everywhere submanifolds of C" with certain 
singularities. Real algebraic geometry studies the set of solutions to systems over real closed fields 
and the solution sets define so-called semialgebraic sets [5] . For an exposition of solution methods 
for polynomial systems, we refer the interested reader to [SS] and [TS]. Note that in our case both 
the problem and its solution set are parametrized by the singular values (Ti > (T 2 > (T 3 > 0 of 
the deformation gradient F € GL''"(3) encoded by the diagonal matrix D — diag((Ti, (T 2 , 0 - 3 ). The 
study of parametrized polynomial systems is an active research area in computational algebraic 
geometry, see, e.g., [S3] and references thereinj^ 

We briefly introduce the Euler-Lagrange equations obtained by taking variations on the matrix 
group SO(3); cf. jSS] p. 28] for details. Let ^ = RA G Tr S0(3) = R ■ so(3) be a direction in the 
tangent space at i? G SO (3). The corresponding directional derivative of the Cosserat shear-stretch 
energy is then 

DrW^^^,^{R-,F).^ = 2y.{sjm{R^F - 1), sym(^^F)) -|-2/ic ( skew(i?^F), skew(^'^E)) 

= (2/r sym([/ — 1) -|- 2^c skew(C7), A'^U) . 

Equating this derivative with zero and noting, as usual, that this equality must hold for all in¬ 
finitesimal rotations A G so (3), we obtain the Euler-Lagrange equations in matrix representation. 
In particular, any critical U '■= R"^F must satisfy 

skew (^{y — y,c) = 0 . (2-17) 


Clearly, the polar factor Rp solves the Euler-Lagrange equations as it symmetrizes U. Thus, Rp 
is always a critical point, see, e.g., m , or ESj. Under certain conditions on F, however, there 
may be non-classical critical points and even minimizers for which U is no longer symmetric! This 
observation lies at the heart of the first collaboration of the present authors [301 ESI and we shall 
meet this phenomenon again in the following; cf. also [80] . 


We have compiled the solution set for the Euler-Lagrange equations in quaternion representa¬ 
tion (2.161 which we have obtained by using Mathematica in Appendix [A] This permits us to 
present the energy-minimizing relative rotations which solve Problem |2 . 5] without further ado. 


Computer Assisted Result 2.6 (Energy-minimizing quaternions for (/r,/ic) = (IjO)). Let 
D = diag(CTi, ( 72 , cr 3 ) with ui > (T 2 > (73 > 0. Then the quaternion representation of the energy- 
minimizing relative rotations for Wi^o{q; D) are given by the following critical points (listed in 
Appendix\^: 


QiAD) 

9n,i(^) 


I3 

±arccos(^^), (0, 0, 1) 


if S 12 := cTi -I- (72 < 2 , 
if S 12 := (7i -I- (72 > 2 . 


(2.18) 


®The present authors are not specialists in (computational) algebraic geometry. Our goal here is to point out 
some interesting references and developments that might be useful for the solution of polynomial systems arising 
also in other applications. 
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Validation. At present, we cannot give a full proof for this result. However, we consider our 
numerical validation to be quite thorough. For an exposition of our analysis of the critical points 
compiled in Appendixj^and the numerical validation of the presented energy-minimizing solutions 
based on extensive random sampling of SO(3) we refer our reader to Section]^ 

One of the main gaps towards a full proof is the question whether the set of critical points computed 
by Mathematica is complete. Note that our extensive validation based on random rotations, which 
exceeds what we can present in a paper by far, does not hint at the existence of additional critical 
points. Solving algebraic problems is the domain where CAS tools such as Mathematica do shine 
brightly. 


Corollary 2.7 (Energy-minimizing relative rotations for (/r,/Xc) = (1,0)). The solutions to Prob¬ 
lem [H are given by the energy-minimizing relative rotations 



Here, the optimal relative rotation angles are given by 


(2.19) 


^Uf) 


0 

±arccos(;;^) 


, if Si2 
, if Si2 


CTl -\- (72 < 2 , 
ai-\- a2>2 . 


( 2 . 20 ) 


In particular, for cti -|- CT 2 < 2, we obtain Rf^^F) = 1. 


The interpretation of the optimal relative Cosserat rotations is the main subject of the next sec¬ 
tion, but in anticipation of this subsequent discussion we remark that the condition cti -I- CT 2 < 2 
characterizes a generalized compressive regime. 


3 Optimal Cosserat rotations, maximal mean planar strain 
and the reduced energy 


All proper rotations of euclidean three-space act in a plane perpendicular to the axis of rotation. 
From this, a continuum model with rotational degrees of freedom inherits a certain planar character. 
In our context, it seems natural to introduce 

Definition 3.1 (Maximal mean planar stretch and strain). Let F G GL~^(n), n > 2, with singular 
values fTi > (72 CTn > 0. We introduce the maximal mean planar stretch u™“p and the 

maximal mean planar strain s™™p as follows: 


ii“”P(E) 

s“”P(E) 


(7i -\- (7-j ai -\- a 2 , 

max-- —- = --- , and 

2 2 


{<7^ 

max — 
i=f3 


1) + (p'i 
2 


1 ) 


= u““P(F) - 1 . 


(3.1) 


Definition 3.2 (Classical and non-classical domain). To any pair of material parameters {p, pc) 
in the non-classical range p, > pc P we associate the following classical domain and non- 
classical domain for the parameter F G GL~^(n) 


Dom[(_^^ := {F G GL+(n) | s™“P(F^,^J < 0} , and 
Dom;:^^ := {F G GL+(n) | s““P(F^,^J > 0} , 

respectively. 


It is straight-forward to derive the following alternative characterizations 

= {Fg GL+(n) | u””P(F) < = {F G GL+(n) | ai -t aa < := , 

l-f f-fc 

= {Fg GL+(n) | u””P(F) > A^.^J = {F G GL+(n) \ ai + a 2 > p^.^, := . 

p - Pc 
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Note that the intersection Dom®,, 
However, the minimizers rpolar^ 


n = {F & GL+(n) | i® empty, 

coincide with the polar factor Rp(F) on this intersection. 
This can be seen from the form of the optimal relative rotations in Corollary |2.7[ In particular, 
for dimension n = 3, we rediscover the following important characterizations of these domains for 
the non-classical limit case {fi, Hc) = ( 1 , 0 ); cf. (2.18): 


Domf g := {F G GL"'‘(3) | Si 2 := + 0-2 < 2} , and 

Dom^^ := {F G GL+(3) | S 12 := (Ti + <12 > 2} . 


Previously, in our Corollary 2.7 we have determined the energy-minimizing relative rotations 


Ri q{D) := argmin , Z?) := argmin sjm{RD) — 1 

ReSO{3) flGSO(3) 


(3.5) 


Let us briefly summarize: for m“™p(F) < 1, i.e., when F G Dom)J^g, we have R^q(D) = 1 which 
corresponds uniquely to the polar factor Rp . The minimizers rpolar)*^g(F) deviate strictly from 
Rp(F) for F G DomJ^g \ Dom^g and are hence non-classical. Further, expressed in terms of the 
maximal mean planar stretch u™™p(F), we obtain the alternative representation 


FifiiF) — 




V 


Ta/1 - 


±Jl- 


0 


Ummp(i?)2 
1 0 


o\ 


u^"'p{F) 

0 


1 / 


(3.6) 


Towards a geometric interpretation of the energy¬ 
minimizing Cosserat rotations rpolar)*^Q(F) in the non- 
classical limit case {fi, fJ-c) = (1,0), we reconsider the spec¬ 
tral decomposition of 17 = QDQ^ from the principal axis 
transformation in Section [T] Let us denote the columns 
of Q € SO(3) by qi G S^, i = 1,2,3. Then qi and 92 
are orthonormal eigenvectors of U which correspond to the 
largest two singular values cti and tT 2 of F G GL“'"(3). More 
generally, we introduce the following 



Definition 3.3 (Plane of maximal strain). The plane of 
maximal strain is the linear subspace 

P'"P(F) := span({(7i,g2}) CK" 

spanned by the two maximal eigenvectors qi,q 2 of U, i.e., 
the eigenvectors associated to the two largest singular val¬ 
ues ui > U 2 >...> cfn of the deformation gradient 
F G GL+(n), n > 2. 

We recall that, due to the parameter reduction [5^1 
Lem. 2.2], it is always possible to recover the optimal rota¬ 
tions for general non-classical parameters p. > Pc ^ 0 

rpolar^ (F) := arg min (F; F) . (3.7) 

flG SO(3) 

from the non-classical limit case {p,Pc) = (1,0)- However, 
we defer the explicit procedure for a bit since it is quite 
instructive to interpret this distinguished non-classical limit 
case first. 


Figure 3.1: A stretch ellipsoid corre¬ 
sponding to (cri,cr 2 ,cr 3 ) = (4, 2, 1/2). 
The plane of maximal strain P™p(F) 
is depicted in blue. The cylinder per¬ 
pendicular to this plane marks the axis 
of rotation 93 T P“p(F) of rpolar=^ (F) 
which corresponds to the eigenvector 
associated with the smallest singular 
value (T 3 = 1/2. The thin blue cylin¬ 
der which bisects the angle enclosed 
by the opening of the ellipsoid corre¬ 
sponds to the polar factor Rp . Each 
of the two outer red cylinders cor¬ 
responds to a non-classical minimizer 
rpolar]*]Q(F). The angle enclosed is the 
optimal relative rotation angle = 
± arccos( ^^^^^ ). This is the major 
symmetry of the non-classical mini¬ 
mizers. 


Remark 3.4 (rpolarfQ(F) in the classical domain). For s“™p(F) < 0 the maximal mean planar 
stretch is non-expansive. By definition, we have F G Dom)]^Q in the classical domain, for which 
the energy-minimizing relative rotation is given by Ri^o{F) = 1 and there is no deviation from the 
polar factor. In short rpolar]*;g(F) = Rp(F). 
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Let us now turn to the more interesting non-classical case F G DomJ^Q. 

Remark 3.5 (rpolari’^g(F) in the non-classical domain). If F G DomJ^g, then by definition 
s'nmp(^) > 0 and the maximal mean planar strain is expansive. The deviation of the non-elassical 
energy-minimizing rotations rpolar)*^g(_F') from the polar factor Rp is measured by a rotation in 
the plane of maximal strain P“p(F) given by Rp(F)^rpolarfg(F) = Q{F)Rf Q{F)Q{F)'^. The 
rotation axis is the eigenvector associated with the smallest singular value 173 > 0 of F and the 
relative rotation angle is given by Pf^iF) = =F arccos (l/u““P(i^)). The rotation angles increase 
monotonically towards the asymptotic limits 


lim 


PtoiF) 


iLtf . 


In axis-angle representation, we obtain 


rUf) 

Rp"^ rpolarf g(F) 


[±arccos(l/u™™P(F)), (0, 0, 1)] , 

[Tarccos(l/M™P(F)), 93] • 


and 


(3.8) 

(3.9) 


Corollary 3.6 (An explicit formula for rpolar)^^^(A')). For the non-elassical limit case (p,, fj-c) = 
(1,0) we have the following formula for the energy-minimizing Cosserat rotations: 


rpolar)^o(F) 


Rp(F) ,if FeDom)J^_o, 

Rp(F)Q(F)i?To(F)Q(F)^ ,if F e Domf’^ . 


For general values of the weights in the non-classical range fi > fXc ^ 0, we obtain 

rpolar;^,^^(^') := rpolarf_o(p;.,^J , (3.11) 

where := F is obtained by resealing the deformation gradient with the inverse of the 

induced scaling parameter := jpzfp- > 0. 


Proof. This is a straightforward application of our equation ( |2.4[ ) which translates relative to 
absolute rotations derived in Section to the optimal relative rotations described in Corollary 
2.7 The second part is non-trivial and follows from [5S1 Lem. 2.2]. ■ 


Note that the previous definition is relative to a fixed choice of the orthonormal factor Q{F) G 
SO(3) in the spectral decomposition of 17 = QDQ^. Further, right from their variational char¬ 
acterization, one easily deduces that the energy-minimizing rotations satisfy rpolar))^^^((5F) = 
Q rpolar)^ (F), for any Q G SO(3), i.e., they are objective functions; cf.Remark 


3.7 


3.6 


indicate a certain tension- 


The domains of the piecewise definition of rpolar]*^g(F) in Corollary 
compression asymmetry in the material model characterized by the Cosserat shear-stretch energy 
Wip{R-, F); cf. Remark 3.15 We can also make a second important observation. To this end. 


consider a smooth curve F{t) : {—e,e) -G GL''"(3). If the eigenvector ^ 3 ( 1 ) G associated with 
the smallest singular value 173 ( 1 ) changes its orientation along this curve, then the rotation axis 
of rpolar]^g(F) flips as well. Effectively, the sign of the relative rotation angle PfgiF) is negated 
which may lead to jumps. This can happen, e.g., if F(l) passes through a deformation gradient 
with a non-simple singular value, but it may also depend on details of the specific algorithm used 
for the computation of the eigenbasis. 


For the classical range /Xc > /r > 0, the polar factor and the relaxed polar factor(s) coincide 
and trivially share all properties. This is no longer true for the non-classical parameter range 
Me ^ M > 0 s-iid we compare the properties for that range in our next remark. More precisely, 
we present a detailed comparison of the well-known features of the polar factor Rp which are of 
fundamental importance in the context of mechanics. 
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Rp(f) 



Figure 3.2: Pitchfork bifurcation diagram for rpolar^^^(F") for fJ. > Hc ^ 0. Let us express the 
energy-minimizers rpolarj^^^(F") in terms of the maximal mean planar stretch of the 

rescaled deformation gradient := For values F G Dom^_^^, we have 0 < 

and the polar factor Rp(F) is uniquely energy-minimizing. In contrast, for F G DomJ^^^, A^_^^ < 
u™mp < ooj there are two non-classical minimizers rpolar*^^(F"). In this regime, the polar factor 
is no longer optimal but it is still a critical point. At the branching point = A^^^^ the 

minimizers all coincide: rpolar“^^(F) = Rp(F’) = rpolar^^^(F'). For /i^ —?► /i, the branching point 
escapes to infinity which asymptotically recovers the behavior in the classical parameter range 

Hc> fJ'> 0. 


Remark 3.7 (Rp(F) vs. rpolar(F) for the non-classical range fJ. > /ic > 0). Let n > 2 and 
F G GL“''(n). The polar factor Rp{F) G SO(n) obtained from the polar decomposition F = Rp(F)C/ 
is always unique and satisfies: 


(Objectivity) Rp( Q • F") 

(Isotropy) Rp (F • Q) 

(Scaling invariance) Rp( A • F) 

(Inversion symmetry) Rp(F“^) 


Q ■ Rp{F) 

Rp(F) • Q 

MF) 

Rp(F)-i . 


(VQ G SO(n)) , 
(VQ G SO(n)) , 
(VA > 0) , 


(3.12) 


The relaxed polar factor(s) rpolar^_^^(F) C SO(n) is in general multi-valued and, due to its vari¬ 
ational characterization, satisfies: 


(Objectivity) rpolar^^^^ {Q-F) = Q- rpolar^_^^(F) (VQ G SO(n)) , 

(Isotropy) rpolar^^^^(F • Q) = rpolar^_^^(F) • Q (VQ G SO(n)) . 

For the particular dimensions k = 2,3, our explicit formulae imply (cf. also Part I \29)) that there 
exist particular instances A* > 0 and F* G GL~''(fc) for which we have 


(Broken scaling invariance) rpolar)^_^^(A* • F*) 7 ^ rpolar(F*) 

(Broken inversion symmetry) rpolar)^^^^(F*“^) 7 ^ rpolar(F*)“^ 


and 


(3.14) 


This can be directly inferred from the partitioning of GL^(fc) = Dom^ 
respective piecewise definition of the relaxed polar factor(s), see Corollary 


,, U Dom))^'f and the 


3.6 


We interpret these broken symmetries as a (generalized) tension-compression asymmetry. 


13 










3.1 The reduced Cosserat shear—stretch energy 


We now introduce the notion of a reduced energy which is realized by the energy-minimizing 
rotations rpolar^ see also Remark 

Definition 3.8 (Reduced Cosserat shear-stretch energy). The reduced Cosserat sheai—stretch 
energy is defined as 


5.1 


W, 


red 


: GL+(n) 






(3.15) 


Besides the previous definition, we also have the following equivalent means for the explicit com¬ 
putation of the reduced energy 


= TTM.Mc(rpolar±^„(F);F), and 

= K:LiD) :=^min D) . 

iiGSO(n) 


(3.16) 


We now approach the computation of the explicit representation of W^®))^(F) by means of the 
equivalent expression For the sake of brevity, we set c = = 1 /u™“p(F) and 

s = \/l — c^. This allows us to write the optimal relative Cosserat rotations in a simple form in 
the computation of 

/ c =Fs 0 \ /cti 0 0 \ / (Ti • c =Fcr 2 ■ s 0 \ 

RfoD = I ±s c 0 j I 0 (72 1 ~ 1 ^ 0-2 • c 0 1 ■ (3-17) 

V 0 0 1/ yo 0 ( 73 / V 0 0 (73/ 


From this, we compute the following symmetric and skew-symmetric parts: 

\ Ai-c -1 

di2 , „ 

2 * 

0 \ 

sym = dm . g 

CT 2 • C — 1 

0 , 

V 0 

0 

(73 - l/ 

/ 0 



skew ( RTqD — 1 ) = ± ^ ■ s 

0 

0 . 

^ V 0 

0 

0 / 


(3.18) 


Lemma 3.9 (The reduced Cosserat shear-stretch energy W™o(l") in terms of singular values). 
Let F G CL“''(3) and (7i > (72 > (73 > 0 the ordered singular values of F. Then the reduced Cosserat 
shear-stretch energy W\°q{F) admits the following piecewise representation 


^TpiF) 


((7i - 1)2 -I- (ct2 - 1)2^ -I- ((73 - 1)2 = lie - 1||^ , if (7i -I- (72 < 2, i.e., F G Domf g , 
i ((7i - (72)2 -h ((73 - 1)2 , */ 0-1 -t (72 > 2, be., F G DomJl^y . 


Proof. The classical piece of the energy is easily obtained by inserting the polar factor Rp(e) into 
the energy. To compute the non-classical piece, we first recall that 


sym(rpolar±(F)^F-l)f = W^®(((F) = Wlf{D) 


sym{R^D — 1 ) 


0 - 1 - 1 - 0-2 


and s = -s/l — ( 


We compute the expression on the right hand side. To this end, we set c = 
again and compute the Frobenius matrix norm of sym(i?^£) — 1) which we have derived in (3.18). 
This gives us 


sym 


(^R^D - 1 


((7iC - 1)^ -I- ((72C - 1)^ + ^ ((7l - (72)^(1 “ C^) -|- ((73 - 1)^ 

^ (4 + ((7i - CT2)^ - 4 c((7i -k (72) + c2((7i + 02^) + ((73 - 1)^ 

^((7i - ( 72 )^ -k ((73 - 1)^ . ■ 


Our next step is to reveal the form of the reduced energy for the entire non-classical parameter 
range > /Tc > 0 which involves the parameter reduction lemma, but we have to be a bit careful. 
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Remark 3.10 (Reduced energies and the parameter reduction lemma). The parameter reduction 
procedure described in JSdl Lem. 2.2] is the key to the minimizers for general non-classical material 
parameters ^ > /ic > 0. It might he tempting, but we have to stress that the general form of the 
reduced energy cannot be obtained by rescaling the singular values Oi i—)■ in the singular 

value representation of . 

Theorem 3.11 as a function of the singular values). Let F G GL''’(n) and ai > a 2 > cr^ > 

0, the ordered singular values of F and let /r > > 0, i.e., a non-classical parameter set. Then 

the reduced Cosserat shear-stretch energy : GL'''(3) —)■ K)}" admits the following explicit 

representation 




p ((cti - 1)2 + ( 0-2 - 1)2 + ((73 - 1 ) 2 ) = p\\U-l\f 

f (cti — (72)2 + /r ((73 — 1)2 + ^ (((7i + (72) — 


2 


’ P 


2 


, F G Dom 
, F G Dom 


c 

NC 


Proof. In order to obtain the classical part of the energy it suffices to insert Rp into the energy. 
For the non-classical piece, we insert the optimal relative rotations Rt u This 


amounts to replace c ^ c = and si->-s = -\/l — c2in our preparatory calculation (3.18). 

This yields the following contributions: 




sym.{R^ 


D-r 


skew(.R±^^i:) - 1) 


^ ^32+M (o’s — 1)^ + ^ (p/i. 


Pc 


Pc 


’12 


’12 


2 


Pv,Uc 


(3.19) 

(3.20) 


Finally, adding only the constant part of the symmetric contribution to the complete contribution 
due to the skew-symmetric part, we obtain 


P 

2 


(P/i, /ic 


-2)2 + 


Pc 2 
2 12 


Pc 2 

2 Pu. Vc 



2 P/i, fj,y) 


Pc 

y 


(^12 - 


2 


pc 2 
^ Pu.Uc ■ 


The last step of the preceding proof is interesting in its own right. 

Remark 3.12 (On pc as a penalty weight). Let us consider the contribution of the skew-term to 
given by 

y ((o-i +(T2) -Pm,/xJ^ 

as a penalty term for F G GL~''(3) arising for material parameters in the non-classical parameter 
range p > pc > 0. This leads to a simple but interesting observation for strictly positive pc > 0. 
The minimizers F G GL~''(3) for the penalty term satisfy the bifurcation criterion 


(7l + (72 — P^,3i,, 

for rpolar)^^^^(F). In this case R^^^^ = 1 which implies that — 1 G Sym(3), i.e., it is 

symmetric. Hence, the skew-part vanishes entirely which minimizes the penalty. In numerical 
applications, a rotation field R approximating rpolar^(F) can be expected to be unstable in the 
vicinity of the branching point ai + a 2 ~ Pu.uc Hence, a penalty which explicitly rewards an 
approximation to the bifurcation point seems to be a delicate property. In strong contrast, for the 
case when the Cosserat couple modulus is zero, i.e., pc = 0, the penalty term vanishes entirely. 
This hints at a possibly more favorable qualitative behavior of the model in that case; cf. JWi- 

3.2 Geometric aspects of the reduced Cosserat shear—stretch energy 

We recall that the tangent bundle rSO(n) is isomorphic to the product SO(n) x so(n) as a vector 
bundle. This is commonly referred to as the left trivialization, see, e.g., [23]. With this we can 
comfortably minimize over the tangent bundle in the following lemma which sets the course for 
our next theorem. 

Lemma 3.13. Let F G Then 

inf \\R'^F — 1 — A\y = min || sym(i?^F — 1)||2 =: min WioiR',F). 

R&SO{3) R&SO{3) 1 1 fl:GSO(3) ’ ' ^ 

^4 G 00 (3) 
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Figure 3.3: Energy isosurfaces of WI'^q considered as a function of the unordered singular values 
cri,(J 2 ,<J 3 > 0 of F G GL'''(3). The displayed contour levels are 0.1, 0.4 and 0.8. On the right, we 
have removed a piece from the non-classical cylindrical parts (red) of the energy level 0.8 which 
reveals the spherical shell of the classical part (green). Note that a computation of these level 
surfaces via Monte Carlo minimization yields the same result (but at a much lower resolution). 

Proof. For all R G SO(3), the infimum of ||skew(i?^F — 1) — A||^ over all skew symmetric A is 
obviously attained at ^ = skew(i?^F — 1). Therefore 

inf Wr'^F-I-A 
fiGSO(3) " 

.A G so(3) 


Since SO(3) is compact and Wifi{R; F) is continuous, the infimum is attained by a minimizer. ■ 


= inf inf ( ||sym(F^F-1-yl)||% ||skew(F'^F-1-A) 
flGSO(3) AG6 o(3) L " ^11 II ' 

= inf I ||sym(i?^F — 1 )-|- inf ||skew(i?"^F — 1) — All^ I 
flGSO(3) I II ' 'll AG6 o(3) II ' ' II J 

= inf ||sym(i?^F — 1)11^ . 
fl;GSO(3) II 11 


The preceding lemma leads us to a nice geometric characterization of the reduced Cosserat shear- 
stretch energy which we find quite remarkable. It might even be useful for the case n > 4 although 
this is somewhat far-fetched. 

Corollary 3.14 (Characterization of as a distance). Let n > 2 and consider F G GL~^(n) 

with singular values cti > CT 2 > ... > cr„ > 0 not necessarily distinct. Then the reduced Cosserat 
shear-stretch energy : CL'''(n) -G K)}" admits the following characterization as a distance: 

W^^^iF) = distLiid {F. SO(n) (1 + so(n))) . (3.21) 

Here, disteuciM denotes the euclidean distance function. 


Proof. First note that 


(^) 


min II sym(F^F - l)f 
R G SO(n) 


= inf |1F(F^F-1-A)f 
flGSO(3) 

A G SO(3) 


inf \\R^F-1-Af 
ReSO(3) 

A G so(3) 


The last step ist justified by the orthogonal invariance of the Frobenius norm || • ||. Carrying out 
the multiplications on the right hand side, we are lead to the conclusion 

inf IlF- i?l - RAf = inf ||F- F(1 + A)f =: distLiid(^’, SO(n) (1 +so(n))) . 

R G oO(3) R G oO(3j 

A G 5o(3) a G so(3) 
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Rem ark 3.15 (Zero reduced energy and tension-compression asymmetry). A sharp look at Lemma 
3.9 is sufficient to see that the 0-energy level of WI'^q preci sely corresponds to singular value 
tuples of the form (s, s, 1), s G [l,oo)[^ In our Figure 3.3 tuples of this type (and permuta¬ 
tions thereof) correspond to the axes of the cylindrical sheets of the isosurfaces. Let us now 
consider X = i?(l -|- A), R G SO(3), A G so(3), which has the sguared singular values 

(cr^,(T|,cr|) = -I-II , 1 -I-ll^ll^ , 1^ Clearly, such a matrix X does not generate any re¬ 

duced Cosserat shear-stretch energy at all - in perfect accord with Corollary \3.14\ Geometrically, 
U{X) := V X'^X induces a homogeneous blow-up (i.e., a rescaling of arbitrary positive magnitude) 
of the plane of maximal strain P“p(X) while preserving the distance of any given point to this 
plane. Furthermore, there is no possibility of similar energy savings in the compressive range for 
F G GL~'’(3) where the classical piece ofW[^Q is active. It seems to us that this makes a good case 
for a quite remarkable type of tension-compression asymmetry. 


3.3 Alternative criteria for the existence of non-classical solntions 

For fi > fie > 0, i.e., for strictly positive fic > 0, the singular radius satisfies ■= > 2. 

We now define a quite similar constant, namely 

Cuwo '■= “ Pi,o = > 0 . (3.22) 

M Me 

Furthermore, we define the e-neighborhood of a set X C relative to the euclidean distance 

function as 

Ne[X) := {Y G I disteuciid(r. A") < e} . 



Figure 3.4: Illustration of a euclidean e-neighborhood of SO(3) C 

Lemma 3.16 (Classical SO(3)-neighborhood for pc > 0). Let fi > pc > 0, F G GL''"(3) and 
Cm Mr TTirr > 0- Then we have the following inclusion 

N.^.JSO{3)) C DomC^^. (3.23) 

In other words, for all F G GL“''(3) satisfying disteuciid(F', SO(3)) = \\U — l|j^ < ^Cm.mc' polar 
factor Rp is the unique minimizer of {R; F). 

Proof. Since distgm,j;(j(F', SO(3)) = ||17 — 1||^ = X]i=i(^i ~ 1)^ f’Y Grioli’s theorem [51], we find 

distLiid(^, SO(3)) < ^ Cm.mc ^ 2 ((ui - 1)2 + (a2 - 1)^ + (ag - l)^) < 

=> 2 ((cTi - 1)2 -I- ((72 - 1)2) < . 

^Technically, our derivation of Lemma |3.9| does not extend to the case of multiple singular values, but the 
characterization as a distance function in Corollary |3. 141 does not have this limitation. 
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(3.24) 


Further, 0 < (a — &)^ = + 6^ — 2ab implies 2 (a^ + 5^) > + &^ + 2ab and it follows that 

(ai - 1)2 + (a2 - 1)2 + 2(ai - l)(a2 - 1) < . 

Completing the square and taking square roots on both sides, we find 

((^^ _ 1) + (^2 - 1))" < Cm. ±|(ai-l) + (a2-l)| < Wc • (3.25) 

Inserting Cm,M e := Pm.a^c “ 2, we obtain ((Ji-1) + ((T 2 - 1 ) < Pm,Mc“ 2 . This implies cti+(J 2 < Pm.Mc 
and hence F G Dom(t „ . ■ 


Note that the preceding proof can be quite easily adapted to the planar case n = 2 presented 

in P5] . 

Lemma 3.17. Let F G SL(3), i.e., det[F] = cricr 2<73 = 1, where cti > (72 > 0-3 > 0 are ordered 
singular values of F, not necessarily distinct. Then 

SL(3) C Domf)( , (3.26) 

i.e., F induces a strictly non-classical minimizer. Equivalently, det)!^] = 1 implies the estimate 

( 71 “t” (72 ^ 2 . 


Proof. The inequality for the geometric and arithmetic mean shows that 

O'! + f72 + 0-3 1 

-^- > ((7l (72 0-3) 3 = 1. 

It follows that ( 7 i + (72 > 3 — (73 which implies the claim for CT 3 < 1. 
(71 > (72 > as > 0 , the case (73 > 1 contradicts our assumption det[F] > 1 . 


(3.27) 

Due to the ordering 


Remark 3.18. If we make the stronger assumption ( 7 i > CT 2 > (73 > 0, we obtain a strict inequality 
( 7 i + (72 > 2. In that case, F G Dom)^Q \ Dom^ is strictly non-classical. 

Corollary 3.19. Let p. > 0, F £ SL''"(3) and assume that ai > a 2 > as > 0. Then 

F G DomNC\DomCo, (3.28) 

i.e., the minimizers rpolarC(F) 7 ^ Rp are strictly non-classical. 

Proof. Since Am,o = 1) h follows that Fm,o = F. Further Pm.o = Pi,o- Thus, we are in the hy¬ 
potheses of the preceding Lemma|3.17|for the case where the inequality is strict, see Remark[3.18| ■ 


3.4 Application 


Let us now give a short application to our previous findings. We consider a so-called volumetric- 
isochoric split for the geometrically nonlinear Cosserat shear-stretch energy. Note that this ma¬ 
terial model appears in a variety of contexts, see, e.g., , and, 

recently [3 1511184j . Further, similar expressions for the strain energy have been considered in 
the context of plate and shell theories, see, e.g., [||51[171Eni[I|[I3[I7]. 

Let us introduce the isochoric projection F H. Fjso := £ SL(3) of the deformation gradient 

F G GL“''(3) which can also be applied to U := Ef"F. With this notation, we obtain 


W{U) = p 


sym 


U 


\det[t/] 


■7711/3 


- 1 


skew 


U 


\det[{7] 


■7711/3 


- 1 


“Cosserat shear—stretch energy” 


h{det[U]) 


“volumetric contribution” 


= p 

= M 


sym ( ^ , - 1 I 

sym - 1}||^ + /(c 


+ Pc 


skew 



F 

det[f]^/^ 



skew (i?^Fiso — 1) 11 ^ + h(det[F]) . 


-I- /i(det[i?"^T]) 
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The results of the previous subsections, allow us to determine the optimal Cosserat rotations for the 
split energy IT([/) = W{R^F) = (i?; Tiso) + /i(det[F]). Note first that the additional 

volumetric contribution /i(det[T']) penalizes volume change by a scalar function h : IR“'' —)■ Kj 
which is constant with respect to R G SO(3). Therefore, this formulation still gives rise to the 
same optimal Cosserat rotations 


rpolar^,^^(Tiso) 


argmin 

: ^iso) 

RGSO{3) 


= argmin 
i?GSO(3) 


+ /i(det[F])} . 


We can now make an interesting observation. To this end, let e > 0 and consider diagonal matrices 
of the type 

/ Pii, + 0 \ 

De := 0 1 0 e SL(3) . 

\ 0 0 (P/i.Atc ~ 1 + ^)~ / 

The required ordering tri > CT 2 > 0-3 > 0 follows from := > ^ = 2 and holds for the 

entire non-classical parameter range fi > pc ^ 0- Obviously, we have 

CTi -b 0-2 = Pt,, /.c + ^ > Pm, Mo ^ Dom|^_^^ \ Dom° . 


Hence, the intersection SL(3) C (DomJ^'^^ \ Dom^^^^) yb 0 is never empty since it contains for 
all e > 0. Furthermore, the associated optimal Cosserat rotations are strictly non-classical, i.e., 
rpolar;^.^^(£>e) ^ Rp(i:)^) = 1. 

Hence, in order to assure that there can be no strictly non-classical optimal Cosserat rotations 
(for whatever reason) one has to consider material parameters from the classical parameter range 
Pc > P > 0. In this case the Cosserat couple modulus pc dominates the Lame shear modulus p 
and Grioli’s theorem assures that the polar factor Rp(F) is always uniquely optimal [5^ . 

In the distinguished limit case pc = 0, the volumetric-isochoric split precludes the previously 
observed tension-compression asymmetry. In this particular scenario. Corollary |3.19| shows that 
the optimal rotations are always non-classical. Since no bifurcation of the optimal rotations occurs, 
there can be no qualitatively different energetic response under tension and compression for /ic = 0 ; 
cf. also m for a discussion of other implications of a zero Cosserat couple modulus. 


Last but not least, we want to mention that our proposed explicit formulae for optimal Cosserat 
rotations may also lead to improved stability and performance in full scale 3D nonlinear finite 
element computations for media with rotational microstructure. We expect them to be especially 
useful for the highly interesting and numerically challenging case of a material with small internal 
length scale Lc > 0. If, in addition, the volumetric contribution is independent of t he r otation 

can be 


(see above), then the optimal Cosserat rotations rpolar)f_^^(F) proposed in Corollary 3.6 
expected to be ideal candidates for the initialization of the Newton-iterations for the 
microrotations i? : D C K" —)■ SO(n), n = 2, 3. 


field of 


Dissection of critical point structure and computational 
validation of optimality 


We recall that our primary objective for the present work is to derive a formula (or algorithm) 
which allows to compute the set of optimal Cosserat rotations rpolar^_^^(F) C SO(3), i.e., the 
rotations which minimize the Cosserat shear-stretch energy W^_^_.(i?;F) for given F G GL'''(3) 
and weights (p, pc) in the non-classical parameter range p > pc > 0. In the first two sections of 
this contribution, we have hopefully convinced our avid reader that it suffices to solve Problem 


1.5 in order to determine the optimal Gosserat rotations which then solve our original Problem 
1.1 However, in order to simultaneously cross-validate our theoretical derivation (this includes the 
parameter reduction presented in Part I [5S1 Lem. 2.2, p. 4]), we have based our final validation 
on Problem 0 This bypasses all simplificatio n st eps which we have used in order to derive the 


formula for rpolarjj (F) proposed in Gorollary 3.6 but is costly due to the large parameter space. 
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Figure 4.1: A rhombic dodecahedron placed in the space of unordered singular values (cti, (T 2 ) o'a) € 
of F G GL^(3) gives rise to a beautiful geometric characterization of the classical and non- 
classical domains Domf q and DomJ^Q . Pick a face and displace it in normal direction while scaling 
it by its distance to the origin. This creates a convex cone with the scaled faces as cross-sections 
which intersects the polytope. The part of the cone inside the polytope corresponds to singular 
values (cri,CT 2 ,cr 3 ) G in Domjg. The exterior part corresponds to DomJ^g. On the picked face 
itself Dom^Q n DomJ^g, the two branches coincide. (This is how we discovered the non-classical 
bifurcation behavior of rpolar(*^Q(F) in dimension n = 3.) 


4.1 Interactive analysis of the critical point structure 


The solution of the Euler-Lagrange equations ( 2.16[ ) with the computer algebra pack¬ 
age Mathematica returns the 32 critical points compiled in AppendixjA} Note that Mathematica au¬ 
tomatically verifies that the obtained symbolic expressions are indeed solutions for Problem |2.5| 


These symbolic solutions give rise to 32 critical branches : Dom(g’^*^) —>■ S^, 1 < i < 32. Note 
that, we can discard 16 of the branches right away since they are redundant. This is due to the an¬ 
tipodal identification of quaternions under the covering map tt : —>■ SO(3). The critical branches 

are associated with the lifted Cosserat shear-stretch energy formulation In particu¬ 

lar, they project to relative rotations parametrized by a diagonal matrix D = diag(cri, cr 2 , erg). In 
what follows, we identify the space of unordered singular values (cti, cr 2 , < 73 ) of F G GL“''(3) with 
Further, we take the liberty to identify diagonal matrices D = diag(CTi,(T 2 ,' 73 ) with points in 
and shall even write F G K^. 


This allows us to write Dom(g(*^) C i = 1,..., 16, for the maximal domain of definition of the 
z-th critical branch]^ If the solution set is complete (cf. Appendix [A| for a discussion), then, up to 
a set of measure zero, we must hav^ 


y Dom(g(*)) = . 

l<i<16 


Initially, still stumbling in the dark, we compared the critical branches by comparing the different 
realized energy levels given by ITi,o(g^*^ ;F), i = 1,..., 16, for random tuples (f7i,f72,CT3) G 
This allowed us to construct a three-dimensional map for the space of unordered singular values 
by associating the index set of the energy-minimizing critical branches at F = (cri, t72, < 73 ) with the 
point F G in the parameter space. We then mapped each of these index sets to a unique color 
and subsequently explored the parameter space visually. This three-dimensional “optimal branch 
map” allowed us to isolate the energy-minimizing critical branches corresponding to rpolarj*^g(F); 
cf. Figure [4T] which describes the structure that appeared. 

Furthermore, we compared the following minimal energy levels 

min IFf o(gW(F);F) , and min Wf g(g; F) , (4.1) 

1<2<16 

using a Monte Carlo approximation of the right hand side (which we describe in detail in the next 
subsection). This allowed us to detect discrepancies which can only arise due to an incomplete 
set of critical branches. Note that during our whole investigation, we never encountered any such 

® Technically, these maximal domains are implicitly defined by the requirement that the critical coefficient tuple 
computed by Mathematica is real-valued at the point D E i.e., (icb), xb), y(^), 2 :(*)) ^ 

^Equivalently, for a given F E GL"^(3) with distinct singular values, at least one of the critical branches gb) ; 
Dom(gb)) i;nust be energy-minimizing since the Cosserat shear—stretch energy attains its minimum on SO (3). 
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Figure 4.2: Uniform sampling of the circle by rejection sampling. First, the unit square [—1,1]^ 
is uniformly sampled (here, with 2000 samples). Then all samples p with ||p|j > 1 in the red domain 
are rejected. Finally, the remaining samples in the unit disk are normalized, i.e., p i—>■ S S^. 

This yields a uniform distribution on the boundary circle S^. Although this approach can be 
generalized to higher-dimensional spheres its performance does not scale to high dimensions. 


discrepancy. This is a strong indication that the set of critical points computed by Mathematica is 
in fact complete, as one would expect. 

Our next step is to turn our previously described approach into a more systematic computational 
validation of the optimality of the proposed candidates rpolar(^_^^(U). 


4.2 Validation of optimality by Monte Carlo statistical sampling 


We now describe a serious computational approach for the validation of the optimality of our 
proposed candidate formula rpolar)^^^ (U). This approach relies on a well-known, rather simplistic, 
but highly useful (in low dimensions) method for the generation of uniformly distributed random 
rotations due to [53] . 


In what follows, we let K := [—1,1]"^ C denote a hypercube of sidelength 2 centered about the 
origin and define := {a: S AT | ||a:|| < 1}, i.e., as the closed unit ball in Further, we let 
Xx denote a uniformly distributed random variable with values in K and introduce X^i as the 
restriction of Xx to the unit ball. Then, X^4 := Xx\^i is uniformly distributed. The restriction 
can be defined by rejection sampling, i.e., we reject all realizations in AT \ B *^ which lie outside of 
the ball, but accept the first realization inside of B"^ C AT; see Figure 4.2 for an example in the 
plane. 


Theorem 4.1 (Rejection sampling for S^). The random variable Xga := ||^”^|| obtained by nor¬ 
malization is uniformly distributed on with respect to the Lebesgue measure on the sphere which 
we denote by dVgs. 


Proof. This is a standard method which performs quite well in low space dimensions, see, 
e.g., [301IM] or [35] and references therein. ■ 


We recall that is a Lie group double cover of SO (3). A uniform distribution on a compact 
Lie group is defined in terms of the (normalized) Haar measure of the group, see, e.g., [2 P- 9]. 
Such a measure is invariant with respect to the left (or right) group multiplication and is unique 
up to a constant multiple. For an introduction to the Haar measure on Lie groups, see, e.g., [331 
p. 179-194]. It is well-known that the Lebesgue measure dVgs is a bi-invariant (non-normalized) 
Haar measure for the Lie group of unit quaternions. 

Theorem 4.2 (Uniformly distributed random variables on SO(3)). Let Xga be a uniformly dis¬ 
tributed random variable on and tt : —>■ SO (3) the covering homorphism defined in plot . 
Then the random variable Also( 3 ) := TroAgs is uniformly distributed with respect to the (normalized) 
Haar measure on SO(3). 

Proof. This is well-known, see, e.g., [S3]- B 
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Essentially, the covering homomorphism tt : —>■ SO(3) induces a bi-invariant Riemannian metric 
on SO(3) via the pullback (•, •)gQ( 3 ) := Note that the bi-invariant metric on SO(3) is 

unique up to scalar multiples since the Lie-algebra so(3) is simple}^ With respect to the pullback 
metric, tt is a local isometry. Further, since tt is a covering map, the pullback of the invariant 
surface volume measure on given by (7r“^)*dVs3 induces an invariant measure, i.e., a Haar 
measure, on SO(3). 

On a sidenote, the use of Riemannian metrics and geodesics on (matrix) Lie groups in applications 
is currently an active research area, since the computational costs of geometric methods are no 
longer prohibitive. For some interesting recent applications to strain measures in mechanics, see, 
e.g., [64] . Another interesting recent usecase for geodesics on the group of unit quaternions is the 
simulation of eye movements, see [72] . 

We now briefly describe the sampling strategies for the computational validation. 


Remark 4.3 (Sampling the group of rotations SO(3)). Based on the method described in Theorem 


4-2. we have computed a set of samples Q C consisting o/4.629.171 uniformly distributed unit 


quaternions. 


Remark 4.4 (Sampling strategy for F S GL^(3)). We have generated a stream of matrices with 
uniformly distributed coefficients Fij G 1 < *, j < 3 and discarded all samples with 

det[F] < Oru From the remaining samples, we have selected the first 1.000 samples in Dom®^^ 
and Dom^_^^, respectively, and collected them in two sets and 


Remark 4.5 (Limitations of the sampling strategy for F G GL^(3)). For performance reasons 
our sampling strategy takes our expectations into account right from the onset. This can be seen as 
a limitation. Further, our validation is inherently limited to compact subsets of GL’*' (3). However, 
this particular strategy, heuristically produces a reasonable resolution for parameters F G GL'''(3) 
in the vicinity of the branching condition ai + a 2 = (af our Figure 4-4'l- Based on the 

predictions of the analysis of our proposed optimal Cosserat rotations presented in Section^ this 
is without doubt the most interesting parameter sector. 


We are now finally in the position to expose our computati onal validation strate gy f or the global 


optimality of the formula rpolar* (F) stated in Gorollary 


3.6 


cf. also Remark 


3.4 and Remark 
rotations. It is 


[3.51 for a short review of the geometric interpretation of the optimal Gosserat 
important to note that the presented validation scheme is based on the lift 

:=W^,^,(7r(<7);F). 

This formulation is based on the original Gosserat-shear stretch energy W^^^^(i?;F), precisely as 
it appears in the statement of Problem JTt} Glearl y, thi s allows to validate the consistency of the 


simplifications leading us to Problem 1_^ in Section l|^^ This approach also implies that the image 


SO(3) corresponds to an absolute rotation Tr{q) = R. 


of the covering homomorphism tt : 

Let us now present our 
Computational Validation 4.6. Let the sample sets Q C and Jy 


-4 


:= 




GL~^(3) be as previously defined and set the numerical tolerance tol = 10 Then for all > 0 
and 0 (which we have tested) the following relation holds: 


VF G : TT ("arg min Wj {q ; F)^ =toi rpolar^ (F) , 

V 9eC / 

± I 


(4.2) 


where Ri =toi Ff min± ||Fi — F*|| < tol, Fi,Ff G SO(3). 


The following procedure is equivalent, but more explicit. It also corresponds more closely to our 
actual implementation: 


VF G : Rp{Ff ■ TT faigminWl {q;F)\ 

\ qeQ J 

VF G F^°^ : Rp(F)'^ - TT ("argmin IT /^^^((7 ; F)^ 

\ gee / 


=toi {1} , 

^toi [/3^^.(F),g3(F)] . 


(4.3) 


^®Note that the Lie algebra 5o(n) is not simple for the exceptional dimensions n = 2,4. 

^^We have also discarded matrices with non-simple singular values, but since these form a set of measure zero 
this case did never arise, as expected. 

^^Note that this also extends to our use of the parameter reduction [H Lem. 2.2, p. 4]. 
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n 



Parameters: 

= ( 1 , 1 ) 

Pi,i = oo 


Figure 4.3: Optimal relative rotation angle obtained from stochastic (Monte Carlo) min¬ 
imization for the classical limit case fi = pc = ^- We observe that the relative rotation angle 
vanishes up to numerical accuracy, since the polar factor Rp(F’) is always optimal in perfect ac¬ 
cordance with Grioli’s theorem, see [69] and |29l Cor. 2.4, p. 5]. More precisely, this corresponds 
to the prediction /3)*^p(cri + a 2 ) = 0. For multiple examples from the non-classical parameter range 
P > Pc > 0, see Figure 14.4 


on page 
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In order to clarify the meaning of the notation =toi, let 
^1 =toi Ri =toi Rp, cf. also Remark 


3.5 


iF),q,iF)] = 


G SO(3), then 


Remark 4.7 (Additional verification by a Riemannian Newton-scheme (W. Miiller)). For some 
selected values of F £ GL^(3), W. Muller (then at Karlsruhe Institute of Technology 155^). has 
verified that the proposed formula rpolar)^^^ (F") is a critical point for the Cosserat shear-stretch 
energy. He successfully approximated our proposed optimal Cosserat rotations up to machine accu¬ 
racy by using a Riemannian Newton-scheme for the solution of the Euler-Lagrange equations on 
SO(3). Perturbations of the starting values of the Newton-iteration did not indicate the existence 
of alternative solutions realizing lower energy levels. 


In Figure 4.4 we present multiple plots of the energy-minimizing relative rotation angles 


obtained by stochastic (Monte Carlo) minimization. We show plots for different values of 
A corresponding, in itself rather uninteresting, plot for the classical limit case {p,Pc) = (1;0) is 
depicted in Figure for direct comparison. Both figures match our expectations raised by Figure 
|3.2 1 very well and the resolution does improve with higher sample counts. It is instructive to com¬ 
pare these figures with the optimal relative rotation angles for optimal planar Cosserat rotations 
presented in Part I of the present contribution, see |5S]. 


5 Conclusion 

The reduced Cosserat shear-stretch energy for which we have finally obtained an explicit 

form in Theorem |3.1I| admits an interesting abstract interpretation in mechanics. In order to 
reveal this, let us first assume that the microrotations R are spatially decoupled. This is the case 
when the length scale parameter Lc in the full Cosserat model, i.e., including a curvature energy 
contribution, is extremely small or zero. Let us furthermore assume that det[F"] = 1, i.e., that the 
amount of volume distortion is negligible and that a specimen H of this material is subjected to a 
given deformation ip : Q ^ f(^) with deformation gradient F := Vp G CL"''(3). Then the total 
reduced Cosserat shear-stretch energy obtained by integration of the local density given by 

f W;)®^^^(A)dV := f mm W^.^„(R;F)dV (5.1) 
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Parameters: 

= ( 1 , 0 ) 

Pi,o = 2 


Parameters: 

Me) = (1, 4) 

Pl:\ = I 


Parameters: 

(m, Me) = (1, 2 ) 

Pi,i =4 


Figure 4.4: Optimal relative rotation angles for multiple non-classical values /r > /ic > 0. 
The angles are obtained by stochastic (Monte Carlo) minimization of W^^^^{R\F). The dashed 
blue curve shows the predicted value for + (T 2 ) and the dashed red line marks the expected 

bifurcation point at For a direct comparison, we provide Figure 4.3 on page 23 which shows 

the classical limit case {p, Pc) = (1,1); see also Figure 3.2 on page 13 for an illustration and a more 


precise description of the bifurcation behavior predicted by our proposed formula rpolar)^ (F). 
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corresponds precisely to the total energetic response which is generated if the field of microrotations 
R in the specimen instantaneously aligns itself with the field of locally optimal Cosserat rotations 
rpolar^ {Vip). It is important to observe that the the field of optimal Cosserat rotations is purely 
induced by the deformation mapping ip on which it depends by loeal energy minimization and does 
not otherwise depend on boundary conditions, exterior forces, etc. 

A Cosserat material which conforms to the previous description can be nicely embedded into a 
classical framework due to G. Capriz, the description of which is one of many shimmering pearls 
to be found in the impressive body of his work on micropolar materials, see, e.g., mnn], and it 
is with delight that we summarize it in a brief 

p. 49 ], 


Remark 5.1 (Continua with latent microstructure in the sense of Capriz). In his paper 
G. Capriz introduces the notion of a continuum with latent microstructure as follows: 


“I say that the microstructure is latent when, though its effects are felt in the balance 
equations, all relevant quantities can he expressed in terms of geometric and kinematic 
quantities pertaining to apparent placements. ” 

Capriz then gives a more precise definition of the properties a latent micro structure needs to satisfy. 
We shall only repeat the first two: “There is no inertia connected with the micro structure. ”, and, 
“There are no exterior body actions on the microstructure. ” In other words, a latent micro structure 
is coupled with a deformation (p in an instantaneous way. 


The reduced Cosserat shear-stretch energy can be considered as the energetic answer of a 

medium with a rotational microstructure that instantaneously reorganizes its field of microrotations 
i? : n —>■ SO(3) as an energy-minimizing rpolar)^^^(T')-field. This is an example for a latent 
microstructure in the sense of Capriz. 

From a more general perspective, a Cosserat continuum can also be considered as a special case 
of a so-called micromorphic model, see, e.g., ismss] and m- Let us, as before, set the length 
scale parameter Lc governing the curvature contribution to zero. We then observe that such an 


approach always leads to an algebraic side condition, in our case it is given by the equation (2.17), 
which replaces the partial differential equation for the micro-distortion field. This is another, 
more general, example of a continuum with latent microstructure in the sense of Capriz, compare, 
e.g., El, due to G. Capriz himself and also |51| . 

Note that in |22] and [2T], the authors - who are apparently unaware of this established and rel¬ 
atively straightforward interpretation - have, in our opinion, recently reintroduced the framework 
of materials with latent microstructure due to G. Gapriz for such micromorphic continuum models 
under the new name of a hyperelastic material with “internal balance” and an “internally balanced 
solid”, respectively. 

We now continue our conclusion with some thoughts on possible generalizations of our present 
results. 

Remark 5.2 (On generalizations to higher dimensions n > 4). Our solution approach is quite 
specifically tailored to dimension n = 3 since it relies on the covering o/SO(3) by the unit quater¬ 
nions C H. It seems reasonable to assume that the particularly simple geometry of lies at 
the root of the explicit solvability of the Euler-Lagrange equations. The so-called Sphere Theorem 
states that the only spheres that admit a connected compact Lie group structure are and S^, 
see, e.g., 14^ p.289]. Thus, for n > 3, there is no hope at all to recover the particularly simple 
constellation we have quite successfully exploited here. Still, there is a generalization of the unit 
quaternions, namely the so-called spin groups Spin(n). These groups are two-fold covers o/SO(n) 
and closely related to Clifford algebras, see, e.g., and In principle, such techniques might 
be appropriate for a generalization of our present results to higher dimensions, but they are out of 
reach for us. 


Although our exact solution approach does not generalize to higher dimensions, it seems obvious 
that the reduced Problem |1.5| is a very good starting point for the solution of Problem |1.1| in 
dimensions n > 4. Given this particular form, it seems very likely that the minimizers in higher 
dimensions can also be characterized in terms of the eigenvectors oi U = QDQ^ and the singular 
values (Ji, 1 < i < n, of F’ £ GL~'’(n). Similar to the rather simplistic random sampling strategy 
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we have employed here, it might certainly be worthwhile to carry out an initial investigation based 
on a suitable Monte Carlo random sampling approach which is suitable for higher dimensions, see, 
e.g., |16j . On a related note, we have to dampen expectations regarding extensions to anisotropic 
formulations. These seem to be completely out of reach, since a reduction to a formulation in 
singular values is then impossible, see [SB] and m- 


Another interesting question which is raised by our findings is whether the maximal mean planar 
stretch and strain “measures”, i.e., u“™p(F) and s™“p(F), as defined in Definition 3.1 - which 
appear to be such natural concepts in our particular context - are just artifacts of our derivation. 
The same holds for the plane of maximal strain P™p(F) introduced in Definition 3.3 Are there real- 
world materials or material models which can be precisely or at least approximately characterized 
by, e.g., slip in the plane of maximal strain P™P(i^)? Currently, we are not aware of any such 
materials or models. 


In good hope that the presented mechanisms and computational strategies will be at least helpful 
for the derivation of closed-form solutions for Problem [lT] in dimensions n > 3 and that these will 
match our proposed formula rpolarJ^^(T") presented in Corollary 
present contribution with a last 
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for n = 3, we conclude our 


Remark 5.3 (Final remark). As regards suitable values of the Cosserat couple modulus p,c > 0, 
our development shows clearly that there are ultimately only 3 values of particular interest, namely 


= 0 , He = p , 


and 


He = +00 . 
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A Appendix (list of critical points) 


We now detail our computer assisted strategy for the computation of the critical points for the Lagrange func¬ 
tion Li,o- We recall that the Euler—Lagrange equations simplify considerably if one (or more) of the quaternion 
coefficients w,x,y or z vanishes. This is reflected in the solution set computed by Mathematica. 

We recall our shorthand notation for sums and differences of singular values for parameters F G GL"^ (n) introduced 
in Section [ 2 ] 

Sij := cTj -h (jj and dij := ai — <jj , j = 1, 2, 3 . 


A.l Computation of critical points of the Lagrange function Li o 

In order to solve the Euler-Lagrange equations in quaternion representation l |2.16| l , we have used the Reduce com¬ 
mand in Mathematica. This command returns 130 critical points for the Lagrange function Li,o^] The solution set 
is enlarged by particular solutions satisfying certain algebraic relations among the parameters i = 1,..., 3. In 
order to exclude these cases, we have made the following assumptions 

(ji ^ ^ Sij := (Tj -|- C 7 j 0 , and dij := cjj — crj 7 ^ 0 , hj = 1, 2,3, i ^ j . (1-1) 

These assumptions can be passed to the Reduce command in the form of a so-called AssumptionList. This is a 
standard procedure. With these assumptions Mathematica successfully symbolically reduces the full solution set 
comprised of 130 critical points to 32 critical points. Note that a strict ordering of the parameters ai > cr 2 > (ts > 0 
implies that all of these assumptions in are satisfied. We shall refer to the 32 branches so obtained as the set of 
generic solutions, since they coincide with the output of the Solve command. Due to this procedure, the possibility 
of multiple singular values of F and the degeneracy of det[E] = 0 are explicitly excluded. 

Remark A.l (Completeness of the solution set). We want to stress that we present manually refined results 
obtained via the computational algebra system (CAS) Mathematica. Currently, we cannot prove that our solution 
set is complete, but this is quite probably the case as our extensive validation shows. 


It might be an interesting challenge for an expert in (computational) algebr aic geometry to prove that the presented 
list of generic solutions to the polynomial system l |2.16| | is in fact complete }f‘^\ 

Note that the computations can also be based on the alternative continuation tt' : H ^ SO(3) which produces 
slightly different Euler-Lagrange equations but yields the same full and generic solution sets. 

In order to describe the set of generic solutions, we introduce two auxiliary coefficient functions 


CA{t) := 



and 



( 1 . 2 ) 


It is important to note that the function CB{t) is only real-valued for t G [2, 00 ). The generic solutions are critical 
tuples of complex variables (ty, ic, y, 2 , A) G C^, i.e., the solutions obtained by Mathematica are complex-valued. It 
turns out that the critical tuples are real-valued if and only if the radicands in in are non-negative which allows 
to deduce corresponding domains of definition for the critical branches. 


Remark A.2 (Covering symmetry). For any critical pair {q,X), the pair (— g, A) is also a critical point for the 
Lagrange function Li^q- T^he associated reduced energy levels are identical. 


Due to this antipodal symmetry, it suffices to restrict the following presentation to 16 solutions. For their exposition, 
we have decided to introduce three categories. 


A. 1.1 Type I (of critical points) 

Characterization: Precisely one of the coefficients w,x,y and 2 of q{w, x^y, z) is non-zero. The solutions are 
independent of the parameters ai, i = 1,2, 3. 

These correspond to the following energy-minmizing relative rotations: 


Ql,! 

:= ( 1 , 0 , 0 , 0 ) = 

[ 0 , any ] 

= I 3 

= Ql 

'?I,2 

:= ( 0 , 1 , 0 , 0 ) = 

[^,( 1 , 0 , 0 )] 

= diag( 1 , - 1 , - 1 ) 

= Q 2 

QI,3 

:= ( 0 , 0 , 1 , 0 ) = 

[tt, ( 0 , 1 , 0 )] 

= diag(-l, 1 , - 1 ) 

= Qs 

4iA 

:= ( 0 , 0 , 0 , 1 ) = 

[tt, ( 0 , 0 , 1 )] 

= diag(-l, - 1 , 1 ) 

= Q 4 ■ 


Note that the relative rotation represented by qi^i = I 3 corresponds to the polar factor Rp(E). 

^^Note that the Reduce function in Mathematica is supposedly guaranteed to compute a complete solution set. 
This is not necessarily the case for the Solve-command which generates only generic solutions, see m- 

^"^Note that an attempt to compute either a Grobner basis or a primary decomposition for the Euler-Lagrange 
equations ( |2.16[ l using the CAS Singular m with competent assistance by R. Vollmert and L. Kastner [ 8 fj] 187) 
have failed (to flnish within a day). The parameter-dependent polynomial system might be non-trivial to solve by 
computer algebra. It seems to us that Mathematica automatically carries out the case distinctions w = 0 V w 0, 
X = OVx 7 ^ 0, etc., since a computation of a Grobner basis does not seem to terminate either. These case distinctions 
simplify the Euler-Lagrange equations, but we can just speculate here. 
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For the Lagrange multiplier A, we obtain the associated critical values 

Aj 1 := 0 , Aj 2 ■— ^2 "^23 ) ^ 1,3 •“ *^3 “1“ *^1 “1“ ^ ^31 : and Aj 4 := + cr^ + 4 S 12 . 

The realized energy levels of the lifted energy are given by: 

WloiQhl ; D) = (<71 - 1)2 + (<72 - 1)2 + (<73 - 1)2 , 

(91.2 ; D) = (<71 + 1)2 + (<72 + 1)2 + (,73 - 1)2 , 

W^i%(9I.3 ; D) = (<71 - 1)2 + (<72 + 1)2 + (,73 + 1)2 , 

^^l“, 0 ( 91,4 ; D) = (<71 + 1)2 + (<72 - 1)2 + (,73 + 1)2 . 

The solutions of the first type are globally defined. 


A. 1.2 Type II (of critical points) 


Characterization: Precisely two of the coefficients x,y,z of q{WyX,y, z) vanish. Further, the solution only depends 
on pairwise sums of the singular values cr^, i = 1, 2, 3. 

These correspond to the following energy-minmizing relative rotations: 

9n,i ■= (ca(si 2), 0, 0, ±cs(si2)) = [±arccos(2/si2), (0, 0, 1)] , 

9n,2 ■= ('2 a(' 523), ±cs(s23), 0, 0) = [it arccos(2/s23), (1, 0, 0)] , (1.6) 

9ii,3 ■= (ca(s3i), 0, ±cb(s3i), 0) = [±arccos(2/s3i), (0, 1, 0)] . 


For the Lagrange multiplier A, we obtain the associated critical values 

, >'u 2 -=‘^h(— -. an<i \f := 4 

\ Si2 J \ S23 / 

The realized energy levels of the lifted energy are given by: 

^l,o(9lI,l < -D) = ^(<^1 - + (<73 - 1)^ , 

^l,o(9n,2 ; -D) = i(<72 - <73)2 + (<71 - 1)2 , 

^l,o(9n,3 ; -D) = i(<73 - <7i)2 + (<72 - 1)2 . 

The solutions of the second type are defined for sij G [2, oo), i,j = 1, 2,3. 



(1.7) 


( 1 . 8 ) 


A.1.3 Type III (of critical points) 

Characterization: The coefficient w vanishes together with exactly one of the remaining coefficients x,y,z of 
q{w, X, y, z). Further, the solution only depends on pairwise differences of the singular values ai, 2 = 1,2, 3. 

These correspond to the following energy-minmizing relative rotations: 

^ma ■“ ca(c^ 12 ), ±cs(di 2 ), 0 ) = I^TT, ± 03 (^ 12 ) 2 , O^j , 

9ni,2 ■“ Ca(c^23); ±^ 5 (^ 23 )) = [tT; (o, C^(d23)“2, ±CB(d23)2^j , (1.9) 

9ni,3 ■“ ^^(—^ 31)5 0 , ±£^(—^ 31 )) = I^TT, d 3 i)“ 2 , 0 , ±cs(—d3i)2^j . 

For the Lagrange multiplier A, we obtain the associated critical values 

•“ 40-3(1-h 0-3) + (s23 — 2 )s 23 , 

'^ni,2 4o-i(l-h 0 - 1 ) + (si2 - 2 )si2 , (1-10) 

A^i_3 ;= 40 - 2(1 + 0 - 2 ) + (s3i — 2 )s3i . 

The realized energy levels of the lifted energy are given by 

^i,o(^iii,i i -^) “ 2 ’ 

^i,o(^iii,2 i“ -(0-2 + 0-3)^ + (ti- h 1)^ , (1-11) 

^i,o(^iii,3 ’ “ 2 ' 

The solutions of the third type are defined for dij G [2, oo), i < j, z, j = 1, 2,3. 

Remark A.3 (On second order conditions). We have succeeded to compute the signs of the principal minors of the 
so-called bordered Hessian ll(^^x-^Li^o{q, X; D) evaluated at the previously presented 32 critical points. This allows 
to carry out an analysis of the second order conditions for local constrained extrema based on an analysis of the 
sign changes of the principal minors of the bordered Hessian. The procedure in the constrained case is different but 
similar to the well-known procedure in the unconstrained case im\33- 
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